home *** CD-ROM | disk | FTP | other *** search
/ CD-ROM Today - The Disc! 5 / CD-ROM Today - The Disc (Issue 5)(November 1994).ISO / mac / Mac shareware / Education / RLaB / examples / contour.r < prev    next >
Text File  |  1994-09-21  |  783b  |  50 lines

  1. #
  2. # Simple contour plotting demonstration
  3. #
  4.  
  5. NX = 20;
  6. NY = 20;
  7. xx = zeros (1, NX);
  8. yy = zeros (1, NY);
  9. for (i in 1:NX) { xx[i] = (i - NX/2)/(NX/2); }
  10. for (i in 1:NY) { yy[i] = (i - NY/2)/(NY/2); }
  11.  
  12. zz = zeros (NX, NY);
  13.  
  14. for (i in 1:NX)
  15. {
  16.   for (j in 1:NY)
  17.   {
  18.     r = sqrt (xx[i]^2 + yy[j]^2);
  19.     zz[i;j] = exp (-r * r) * cos (2*pi*r);
  20.   }
  21. }
  22.  
  23. # Sin - Cos surface
  24.  
  25. x1 = -3:3:.2;
  26. y1 = -3:3:.2;
  27. z1 = zeros (x1.n, y1.n);
  28.  
  29. for (i in 1:x1.n)
  30. {
  31.   for(j in 1:y1.n)
  32.   {
  33.     z1[i;j] = sin(y1[j]) * cos(x1[i]);
  34.   }
  35. }
  36.  
  37. # Now create some plots
  38.  
  39. xlabel("X-Axis");
  40. ylabel("Y-Axis");
  41. ptitle ("Contour Demonstration");
  42.  
  43. plmesh (<< x = xx; y = yy; z = zz>>);
  44. pause();
  45. plcont (<< x = xx; y = yy; z = zz>>);
  46. pause();
  47. plmesh (<< x = x1; y = y1; z = z1>>);
  48. pause();
  49. plcont (<< x = x1; y = y1; z = z1>>);
  50.